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1.  Introduction 


Interference  cancellation  is  a  technique  mainly  used  to  reduce  co-channel 
interference.  The  basic  scheme  for  interference  cancellation  is  shown  in  Fig.  1 .  This 
figure  shows  that  a  desired  signal  is  corrupted  by  interference.  The  primary  signal 
contains  both  the  desired  signal  s(n)  and  the  noise  component  n2(n).  The  noise 
component  is  a  function  of  the  reference  passing  through  a  channel  with  an  impulse 
response  w(n).  A  reference  of  the  noise  component  ni(n)  is  then  fed  to  a  finite 
impulse  response  (FIR)  filter  that  adaptively  tunes  using  the  difference  between  the 
adaptive  filter  output  and  the  primary  signal.  The  impulse  response  of  the  adaptive 
filter  will  try  to  match  the  impulse  response  w(n).  If  the  output  of  the  adaptive  filter 
is  equal  to  the  noise  component  of  the  primary  signal,  then  the  error  of  the  adaptive 
filter  is  equal  to  the  desired  signal  s(n). 


L 


Fig.  1  Basic  interference  cancellation  scheme 

2.  Development 


2.1  Hardware  Platform 

To  reduce  development  time,  commercial-off-the-shelf  (COTS)  radios  sold  by 
National  Instruments  subsidiary  Ettus  Research  were  used  as  a  development 
platform.  Interference  cancellers  at  minimum  require  2  inputs  and  1  output.  Ettus 
Research  sells  a  number  of  software-defined  radios  (SDRs)  that  are  capable  of 
meeting  this  requirement.  These  SDRs  are  also  used  with  the  GNU  Radio  software 
development  kit  (SDK),  an  open  source  platform  that  is  available  for  free  in  Linux. 

The  SDR  chosen  for  testing  is  the  universal  software  radio  peripheral  (USRP)  B2 1 0 
and  is  shown  in  Eig.  2. 
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Fig.  2  Ettus  USRP  B210 

According  to  Ettus  Research,  the  B210  has  the  following  specifications: 

1)  Fully  integrated,  2-ehannel  device  with  continuous  radio  frequency  (RF) 
coverage  from  70-6  GHz. 

2)  Full  duplex,  multiple  input  multiple  output  (MIMO)  (2  transmit  and  2 
receive)  operation  with  up  to  56  MHz  of  real-time  bandwidth  (61.44  MS/s 
quadrature) 

3)  FISB  3.0  conneetivity 

4)  GNFl  Radio  and  OpenBTS  support 

5)  Open  and  reconfigurable  Spartan  6  XC6SFX150  field-programmable  gate 
array  (FPGA) 

6)  Analog-to-digital  converter  (ADC)  and  digital-to-analog  converter  (DAC) 
resolution  of  12  bits. 

7)  Receive  noise  figure  of  less  than  8  dB 

Based  on  the  USRP  B210  specifications,  the  theoretical  interference  cancellation 
performance  achievable  is  6  dB  x  12  bits  =  72  dB. 
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2.2  Software  Platform 


Signal  processing  algorithms  were  developed  using  the  GNU  Radio.  GNU  Radio 
is  a  free  and  open  souree  SDK  that  provides  signal  processing  bloeks  to  implement 
software  radios.  It  is  widely  used  in  hobbyist,  academie,  and  eommereial 
environments  to  support  both  wireless  communieations  researeh  and  real-world 
radio  systems.  Using  the  GNU  Radio  simplifies  the  development  process  and 
allows  testing  of  different  of  hardware  platforms  with  only  minor  ehanges  in 
eoding. 

The  signal  proeessing  algorithms  were  developed  in  Ubuntu  14.04  and  the  GNU 
Radio  SDK  version  is  3.7.7. 1.  The  universal  hardware  driver  (UHD)  version  used 
for  the  development  is  003.008.004. 

2.3  Adaptive  Filter  Algorithm 

The  adaptive  algorithm  that  was  used  to  ehange  the  FIR  filter  coeffieients  is  the 
normalized  least  means  squares  (NLMS).  NLMS  is  real-time  recursive  algorithm 
that  adapts  the  FIR  filter  for  every  new  sample  received.  The  algorithm  minimizes 
the  mean  square  error  of  the  difference  between  the  noise  received  by  the  primary 
signal  as  shown  on  Fig.  1  and  the  reference  noise. 

The  pseudoeode  for  the  normalized  least  means  square  algorithm  implemented  is 
shown  below. 

For  eaeh  n 


y(^)  = 

e(n)  =  d(n)  -  y(n)  (2) 

Wj+i(n  +  1)  =  Wj(n)  -  \i  "  (3) 

||x(ti)||  +a 

where  w  is  the  FIR  filter  weight,  x  is  the  noise  referenee  input,  d  is  the  primary 
input  that  eontains  the  desired  signal  as  well  as  the  noise  component,  p,  is  the 
learning  rate,  and  a  is  a  small  number  that  is  greater  than  0.  a  ensures  that  the 
algorithm  is  stable  in  the  event  |  |x(n)  1 1  =0.  The  algorithm  was  designed  to  have 
an  a  of  0.001,  whieh  guarantees  stability  but  also  minimizes  its  error  contribution 
to  the  updated  filter  eoefficient. 

Equation  1  is  the  output  of  the  FIR  filter.  The  seeond  equation  finds  the  differenee 
between  the  primary  input  and  the  output  of  FIR  filter.  If  the  output  of  the  adaptive 
filter  is  equal  to  the  noise  eomponent  within  the  primary  signal,  then  the  error 

Approved  for  public  release;  distribution  unlimited. 


3 


should  only  contain  the  desired  signal.  If  not,  Eq.  3  eonstantly  tunes  the  FIR  filter 
weights  at  a  rate  of  |i.  The  condition  for  |i  must  satisfy  the  following: 

0  <  g  <  1 

The  learning  rate  can  be  arbitrarily  ehosen  and  largely  depends  on  the  type  of 
environment  in  which  the  application  is  being  used  and  what  type  of  signals  are 
involved.  Theoretieally,  small  values  of  p,  inerease  the  amount  of  time  the  algorithm 
takes  to  reach  the  optimal  solution  but  produee  fewer  estimation  errors.  Larger 
values  have  the  opposite  effect  and  increase  the  stability  of  the  algorithm. 

The  algorithm  was  developed  in  GNU  Radio  using  C++.  The  code  is  available  in 
the  Appendix. 

3  Performance 


3.1  Test  Setup 

The  performance  of  the  interference  cancellation  system  was  evaluated  in  a  closed- 
loop  coaxial  cable  environment  for  initial  testing.  The  diagram  of  the  test  setup  is 
shown  in  Fig.  3. 


_ I _ 

Laptop 


Signal  Generator 


Fig.  3  Interference  cancellation  test  setnp 

The  signal  generator  generates  a  single  tone  signal  that  is  fed  into  an  RF  splitter. 
The  outputs  of  the  RF  splitter  is  then  connected  to  each  of  the  receive  inputs  of  the 
USRP  B2I0.  The  USRP  B2I0  feeds  all  the  samples  to  the  laptop  into  the  GNU 
Radio  framework  for  signal  processing.  The  laptop  also  shows  the  performance  of 
the  adaptive  filter  signal  processing  block  written  by  plotting  the  Fourier  transform 
in  real  time.  The  spectrum  plots  were  eaptured  to  measure  cancellation 
performance.  The  interference  canceller  was  tested  for  various  filter  lengths  from 
1,  2,  3,  and  5. 
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The  splitter,  Mini-Cireuits  ZFRSC-42-S+,  and  cable  path  incurs  a  loss  of  7  dB  at 
300  MHz.  When  the  B210  is  configured  in  MIMO  configuration,  the  B210  incurs 
another  8  dB  of  loss. 

3.2  Test  Parameters 

The  following  test  settings  were  used  to  configure  the  signal  generator  and  the 
interference  canceller. 

Signal  Generator: 

1)  Frequency:  300  MHz 

2)  Signal  Type:  Single  Tone 

3)  Modulation  Off 

4)  RF  Amplitude:  -60  dBm 
GNU  Radio  Settings: 

Adaptive  Filter  Parameters 

1)  Sampling  Rate:  64  kHz 

2)  Filter  Length:  1,  2,  3,  and  5 

3)  Learning  Rate  p,  =  0.01 
UHD  Source 

1)  Sampling  Rate:  64  kHz 

2)  Center  Frequency:  300  MHz 

3)  Clock  Frequency:  30.72  MHz 

3.3  Test  Results 

The  amount  of  cancellation  performance  was  recorded  by  measuring  the  spectrum 
at  the  output  of  the  adaptive  filter.  The  measurements  were  taken  at  various  filter 
lengths  of  1,  2,  3,  and  5.  As  shown  in  Figs.  4,  5,  6,  and  7,  the  best  cancellation 
performance  achieved  is  when  the  adaptive  filter  length  is  1.  Figure  4  shows  that 
the  cancellation  performance  is  achieved  is  120  dB.  The  result  is  somewhat 
misleading  because  the  noise  floor  of  the  B210  is  at  -125  dBm,  and  therefore,  the 
realistic  cancellation  performance  is  —60  —  8  —  7  —  (—  125)  =  50  dB.  The 
output  of  the  filter  also  has  a  fiat  spectrum,  indicating  that  the  filter  does  not 
generate  errors  in  performance.  Filter  lengths  2,  3,  and  5  have  noise  suppression 
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Power  (dB)  Power  (dB) 


levels  of  40  dB  at  the  eenter  frequeney  but  generate  an  additional  frequeney 
response  that  inereases  the  noise  in  other  bands.  The  amount  of  noise  generated  ean 
be  as  large  1 5  dB  above  the  noise  floor. 
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Fig.  4  Adaptive  filter  length  of  1 
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After  Cancellation 
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Fig.  5  Adaptive  filter  length  of  2 
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Fig.  6  Adaptive  filter  length  of  3 
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After  Cancellation 
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Fig.  7  Adaptive  filter  length  of  5 

The  results  of  the  performance  of  the  interference  canceller  is  not  surprising.  The 
filter  length  of  1  performs  best  because  the  test  setup  is  a  closed-loop  setup  with 
cables.  The  adaptive  filter  weights  represent  the  impulse  response  of  the  test  setup. 
In  a  closed-loop  setup,  the  signal  does  not  generate  any  delayed  reflections  that 
lengthen  the  impulse  response  of  the  propagation  path.  It  is  also  possible  that  the 
reflections  generated  are  below  the  noise  floor  of  the  radio  and  cannot  be 
represented  by  the  limited  dynamic  range  of  the  radio.  When  the  filter  length  is 
greater  than  1  in  a  closed-loop  system,  the  adaptive  filter  generates  estimation  errors 
that  increase  the  noise  at  the  output. 

4.  Summary  and  Conclusion 

The  performance  of  an  interference  cancellation  prototype  developed  using  GNU 
Radio  and  the  Ettus  radio  B210  has  been  measured.  The  system  is  able  to  achieve 
cancellation  levels  up  to  50  dB  when  the  filter  length  is  properly  chosen. 

The  performance  of  the  system  can  be  further  improved  by  using  more  expensive 
models  from  Ettus  Research.  The  biggest  limitation  on  cancellation  performance 
ADCs  is  the  B210.  Other  Ettus  radios  such  as  the  N210  and  the  X310  have  high 
performance  14-bit  ADCs  and  a  greater  dynamic  range. 
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Appendix.  GNURadio  C++  Code 
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A-1  Adaptive  Filter  C++  Header  File 

#ifndefINCLUDED_ADAPTIVEFILTER_NLMS_IMPL_H 
#define  INCEUDED  ADAPTIVEFILTER  NEMS  IMPE  H 

#include  <AdaptiveFilter/nlms.h> 

#include  <gnuradio/ filter/ fir_filter.h> 

#include  <gnuradio/math.h> 

#include  <vector> 

#include  <stdexcept> 

#include  <gnuradio/gr_complex.h> 

namespace  gr  { 
namespace  AdaptiveFilter  { 

class  nlms  impl :  public  nlms,  filter: :kernel::fir_filter_ccc 

{ 

private: 

gr  complex  d  error; 
gr  complex  d  FilterOutput; 
gr  complex  d  norm; 

std: : vector<gr_complex>  d  new  taps; 
fioat  d_mu; 
fioat  d  FilterEen; 

bool  d  updated; 

protected: 

gr  complex  error(const  gr_complex  &in); 

void  update_tap(gr_complex  &tap,  const  gr  complex  &in); 

public: 

nlms_impl(fioat  mu,  int  FilterEen); 

~nlms_impl(); 

void  set_taps(const  std::vector<gr_complex>  &taps); 
std::vector<gr_complex>  taps()  const; 


{ 

} 


fioat  gainQ  const 
return  d_mu; 


Approved  for  public  release;  distribution  unlimited. 
12 


void  set_gain( float  mu) 

{ 

if(mu<  O.Of  II  mu>  l.Of)  { 

throw  std::out_of_range("NLMS::set_gain:  Gain  value  must  be  in 

[0,1]"); 

} 

d_mu  =  mu; 


//  Where  all  the  action  really  happens 
int  work(int  noutput_items, 

gr_vector_const_void_star  &input_items, 
gr_vector_void_star  &output_items); 


}  //  namespace  AdaptiveFilter 
}  //  namespace  gr 

#endif /*  INCLUDED  ADAPTIVEFILTER  NLMS  IMPL  H  */ 
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A-2  Adaptive  Filter  C++  Source  File 


#ifdef  HAVECONFIGH 
#include  "config.h" 

#endif 

#include  <gnuradio/io_signature.h> 
#include  "nlms_impl.h" 

namespace  gr  { 
namespace  AdaptiveFilter  { 

using  namespace  filter:  :kemel; 

nlms::sptr 

nlms::make(float  mu,  int  FilterFen) 

{ 

return  gnuradio:  :get_initial_sptr 
(new  nlms_impl(mu,  FilterFen)); 

} 


/* 

*  The  private  constructor 

*/ 

nlms_impl::nlms_impl(float  mu, int  FilterFen) 

:  gr::sync_block("nlms", 

gr::io_signature::make(2,  2,  sizeof(gr_complex)), 
gr::io_signature::make(l,  1,  sizeof(gr_complex))), 
fir_fdter_ccc(  1 ,  std:  :vector<gr_complex>(FilterFen,  gr_complex(0,0))), 
d_new_taps(FilterFen,  gr_complex(0,0)), 
d_updated(false),  d_error(gr_complex(0,0)) 

{ 

set_gain(mu); 
if  (FilterFen  >  0){ 

d_taps[0]  =  1.0; 
d_new_taps[0]=  1.0; 

} 

fir_filter_ccc : :  set_taps(d_new_taps); 
set_history(FilterFen); 

} 


/* 

*  Our  virtual  destructor. 

*/ 

nlmsimpl:  :~nlms_impl() 

{ 
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} 

void 

nlms_impl::set_taps(const  std::vector<gr_complex>  &taps) 

{ 

d_new_taps  =  taps; 
dupdated  =  true; 

} 

std: :  vector<gr_complex> 
nlms_impl::taps()  const 

{ 

return  d  taps; 

} 

void 

nlms_impl::update_tap(gr_complex  &tap,  const  gr  complex  &in) 

{ 

tap  =  std::conj‘(tap); 

tap  +=  d_mu*std::conj(d_error)*in/(d_norm+.0001f); 

} 

grcomplex 

nlms_impl;:error(const  gr  complex  &sig) 

{ 

gr  complex  error  =  sig  -  d  FilterOutput; 
return  error; 

} 

int 

nlmsimpl:  :work(int  noutputitems, 

gr_vector_const_void_star  &input_items, 
gr_vector_void_star  &output_items) 

{ 

const  gr  complex  *ref  =  (const  gr  complex  *)  input_items[0]; 

const  gr  complex  *sig  =  (const  gr  complex  *)  input_items[l]; 
gr  complex  *out  =  (gr  complex  *)  output_items[0]; 

if(d_updated)  { 
dtaps  =  dnewtaps; 
set_history(d_taps.size()); 
d_updated  =  false; 

return  0;  //  history  requirements  may  have  changed. 

size_t  k,  1  =  d_taps.size(); 
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} 


size_t  m  =  d_new_taps.size(); 

for  ( int  i  =  0;  i  <  noutput  items;  i++)  { 
dnorm  =  gr_complex(0.0f,0.0f); 
dFilterOutput  =  gr_complex(0.0f,0.0f); 
dFilterOutput  =  fdter(&ref[i]); 

derror  =  error(sig[i]); 

//Calculate  the  norm 
for(  k  =  0;  k  <  1;  k++)  { 

d  norm  +=  ref[i  -  k]*std::conj(ref[i  -  k]); 

} 

derror  =  error(ref[i]); 
out[i]  =  derror; 

for(k  =  0;  k  <  1;  k++)  { 

//  Update  tap  loeally  from  error. 
update_tap(d_taps[k],  ref[i-k]); 

//  Update  aligned  taps  in  fdter  object. 
fir_fdter_ccc::update_tap(std;;conj(d_taps[k]),  k); 

} 

} 

return  noutput_items; 

} 

}  /*  namespaee  AdaptiveFilter  */ 

}  /*  namespaee  gr  */ 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


ADC 

analog-to-digital  converter 

ARL 

US  Army  Researeh  Laboratory 

COTS 

eommereial-off-the-shelf 

DAC 

digital-to-analog  eonverter 

FIR 

finite  impulse  response 

FPGA 

field-programmable  gate  array 

MIMO 

multiple  input  multiple  output 

NLMS 

normalized  least  means  squares 

RF 

radio  frequency 

SDK 

software  development  kit 

SDR 

software-defined  radio 

UHD 

universal  hardware  driver 

USRP 

universal  software  radio  peripheral 
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